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Lewis Research Center 

SUMMARY 

The resonant modes of a dielectric sphere immersed in a warm collisional plasma 
are determined for various electron densities, electron collision frequencies, and elec- 
tron temperatures. It is found that variations in the electron density, collision fre- 
quency, and temperature cause sufficiently large changes in the resonant frequencies and 
damping factors of these modes so that a dielectric sphere may be useful as a plasma 
diagnostic tool. 


INTRODUCTION 

In the many fields where plasmas are used, it has become increasingly important to 
determine accurately the densities, temperatures, and collision rates of the various 
species of particles that constitute the plasma. High-frequency electromagnetic waves, 
or microwaves, have been found to be particularly useful as diagnostic tools for deter- 
mining the properties of the electrons. 

Wait (refs. 1 and 2) has recently pointed out that a spherical region in a plasma that 
is devoid of charged particles may be useful as a plasma diagnostic tool. This spherical 
region acts as an electromagnetic resonator with its resonant frequencies and damping 
factors being functions of the plasma parameters . Such a region can be achieved in 
practice by immersing a dielectric sphere, either solid or hollow, in the plasma. 

This report presents an analysis of the resonant modes in a solid dielectric sphere 
that is immersed in a warm collisional plasma. The resonant frequencies and damping 
factors of several modes are computed for various electron densities, electron temper- 
atures, and electron collision frequencies to determine the feasibility and limitations of 
using this technique as a plasma diagnostic tool. The associated problem of exciting the 
electromagnetic fields in the dielectric sphere is not discussed herein. 



BASIC EQUATIONS 


The dynamics of the plasma will be described by using the linearized hydrodynamic 
and Maxwell equations for a single fluid model of a warm plasma medium (refs. 1 to 3). 
These equations are 
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icop + m n u V - v = 0 
L e o 
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(2) 

Vxh = n qv + iwe e 
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where 

e electric field intensity 

h magnetic field intensity 

v electron fluid velocity 

p perturbation electron pressure 

m mass of electron 

e 

q electronic charge 

u acoustic velocity in electron gas 

n Q equilibrium electron number density 

e Q electric permittivity of free space 

/x magnetic permeability of free space 

co angular frequency 

v effective electron collision frequency for momentum transfer 

It has been assumed that all nonstationary quantities vary with time as e* w * and that c o 
is sufficiently large so that only the electron motion has to be considered. 

Equation (1) is a combination of the continuity and energy conservation equations for 
the electron gas. For the case where the electron gas obeys the equation of state for an 
ideal gas, the acoustic velocity u is given by 


2 



u = 


1/2 


(5) 



where 

k Boltzmann constant 

T g electron temperature 

y ratio of specific heats with constant pressure to constant volume of electron gas 

Equation (2) is a statement of the conservation of momentum for the electrons, and 
equations (3) and (4) are Maxwell equations where the electric current density has been 
set equal to the electron convection current n Q qv. 

The basic equations can be solved more easily if they are expressed in the following 
form: 
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is the effective relative dielectric constant of the plasma 



V* 

e o m 


( 6 ) 

(7) 

( 8 ) 

(9) 


3 


where 


= k 2 e 

op 


k 2 -^ 

V-; 


u 



i 



e 


P 


co 

c 


c velocity of light 

Equations (6) to (9) show that the fields in a warm plasma medium can be determined 
by first solving for the magnetic and pressure fields, h and p, from equations (6) and 
(7), and then computing the electric and velocity fields, e and v, from equations (8) 
and (9). Equations (6) to (9), of course, apply in a dielectric medium if v and p are 
set equal to zero and if is set equal to the relative dielectric constant of the medium. 
Thus, in a dielectric 


V x V xh - Kk 2 h = 0 

e = — - — V xh 

ic % K 

where K is the relative dielectric constant of the medium . 
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MODEL AND BOUNDARY CONDITIONS 


The model to be 
sphere, of radius a 


considered, shown in figure 1, consists of a solid dielectric 
and relative dielectric constant K, that is immersed in an infinite 


Warm homogeneous 
collisfonal plasma 
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^Dielectric sphere 

Figure 1. - Model of dielectric sphere. 



homogeneous warm plasma. It will be assumed that the dielectric material has no elec- 
trical losses. 

Before a solution for the field can proceed, boundary conditions must be specified 
for the field components at the plasma-dielectric interface r = a and at infinity. At the 
plasma-dielectric interface, the components of the electric and magnetic fields that are 
tangent to the surface r = a must be continuous. Thus, 


e x 5 


= e x a. 


r=a 


r=a 


( 12 ) 


h xa r 


_ = h x a 

r=a 


r=a 


(13) 


where a r is the unit vector in the radial direction. A boundary condition for the pres- 
sure or velocity fields must also be specified at r = a to ensure that the fields will be 
unique (ref. 4). The boundary condition that will be used requires the normal component 
of the electron fluid velocity to vanish. Thus, 


v ■ a. 


r=a 


(14) 


This rigidity boundary condition is obviously an approximation. The validity of this con- 
dition has been the subject of much debate (refs. 5 to 7), since the basic equations (eqs. 
(1) to (4)) do not apply in the sheath region that exists at a plasma-dielectric interface 
and since, in many cases, this boundary condition overconstrains the solution in the limit 
of zero electron temperature. A more realistic approach is to approximate the sheath 
region with a thin layer of free space over the surface of the dielectric sphere. To sim- 
plify the analysis, however, this free-space layer is omitted. A discussion of its effect 
on the resonant frequencies is presented in the section NUMERICAL RESULTS. Since 
this analysis is concerned with resonance phenomena localized in the vicinity of the di- 
electric sphere, the proper boundary condition at infinity is to require the field to corre- 
spond to a divergent traveling or decaying wave as r approaches 00 . 

It will be convenient in the following analysis to express all boundary conditions in 
terms of the magnetic and pressure fields. This can be done by combining equations (8), 
(9), and (11) with equations (12) to (14) to give 
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FORMAL SOLUTION FOR RESONANCE EQUATIONS 

The general solution of equations (6) and (10) for the magnetic field can be shown 
(ref. 8) to consist of linear combinations of V x (a r ri//) and V x V x (a r ri//), where the 
generating function ^ satisfies the equations 

V 2 <// + Kk^ = 0 r < a (20) 

vfy + k 2 i// = 0 r > a (21) 

Solutions of equations (20) and (21) that are appropriate for the problem under considera- 
tion are linear combinations of the eigenfunctions 4 / nm where 

Km = e im< ^P™(cos 0)j n (VK k Q r) r < a (22) 

Km = e im< ^P™(cos 0)k n (ik e r) r > a (23) 

The functions P™(cos 0) are the associated Legendre polynomials, and the k Q rj 

and k n (ik g r) are the spherical and modified spherical Bessel functions, respectively. 

The spherical Bessel functions can be constructed by using the following recurrence 
equations: 
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The spherical Bessel functions of the first kind must be used in the solution for the re- 
gion r < a since the field must be bounded at the origin. The modified spherical Bessel 
functions must be used for r > a, where the branch of k g , where Rek g > 0, must be 
selected to ensure that the field corresponds to a divergent wave, in order for the field 
to satisfy the proper boundary condition at infinity. If Rek g = 0, the branch where 
Imk g < 0 should be selected so that the field decays in amplitude as r approaches 00 . 
The general solution for the magnetic field can now be constructed by forming linear 

combinations of the quantities V x (a r xL> ) and V x V x (2 ri 1/1 with the result 
^ ' r nm' r nm' 
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where A B„ . C . and D are constants to be determined, and the prime de- 

nm 7 nm 7 nm 7 nm 7 r 

notes differentiation with respect to the coordinate r or 0. The summations with re- 
spect to m in the previous equations need only be taken from -n to n since the asso- 
ciated Legendre polynomials P™(cos 0) are zero when |m| > n. 

The general solution of equation (7) for the pressure is of the form 

p -t t e im<p p m( cos 6)k n (ikpr) r > a (26) 

n=0 m=-n 

where the P nrn are constants to be determined. The branch of where Rekp>0, 
or Imkp < 0 if Rekp = 0, must be selected to ensure that the pressure satisfies the 
proper boundary condition at infinity. The pressure is, of course, zero in the region 
r < a. 

Relations among the five sets of constants B„ m , C„ rv ,, D„ m , and P„ m can 

° nm’ nm’ nm’ nm’ nm 

be determined by substituting the general solutions for the magnetic and pressure fields 
into the five boundary condition equations, equations (15) to (19). After a considerable 
amount of work, the five sets of equations can be expressed in the following forms: 
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Equations relating the quantities ^/ nm , ^m’ ^nm’ e * c * can be obtained by using 
the orthogonal properties of the exponential functions and the associated Legendre poly- 
nomials. Each term in equations (27) to (31) depends on the coordinate <p only through 
the factor e im< ^; therefore, these equations are also valid with the summation with re- 
spect to m eliminated because the function e im< ^ is orthogonal to if li= m. 

Equations (27) and (28), with the summation over m and the factor e lm< ^ removed, 
can now be combined. Multiplying each term in equation (27) by imP^cos 0)/sin 6 
and each term in equation (28) by P^ m (cos 8) and then adding yield 
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Next, each term in equation (32) is multiplied by sin 8 and integrated with respect to 8 
from 0 to 7 r. Since 
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and 


f* m[p“(cos 0)P™(cos 0)] d6 = 0 


the only way the integrated form of equation (32) can be satisfied is for J nm = # nm 
for all values of n and m except for the special case where n = m = 0. If each term 
in equation (27) were multiplied by Pl m (cos 6) and each term in equation (28) by 


lm 


p“(cos 0) 


sin 6 


the result would be 


nm 


nm 


Thus, equations (27) and (28) are 


satisfied only when ^ nm = ^ nrn and ®' nm = ®nm' Since equations (29) and (30) are of 
the same form as equations (27) and (28), it follows that sd' = <£' - & and 

# nm = ® nm - Equation (31) simply requires that ^ = n(n + l)^ nm - Thus, the de- 

sired relations between the quantities j^ nm , sd' nm , etc. are 


sd = V 

nm nm 
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= @ T 
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rJ T - Cf f _ 

^ nm nm nm 

(35) 
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<m = n < n + 

(37) 


For the special case where n = m = 0, it can be shown that the quantities sd , sd' , 

nm’ nm’ 

etc., are completely arbitrary, except for and &' nm , which must be zero. 

Comparing equations (33) to (37) with the definitions for the quantities ^ nm , ^m’ 

etc. reveals that the equations relating the constants A , B , C , D . and P 

nm’ nm’ nm’ nm’ nm 

divide into two sets with one set involving only B _ and D ^ and the second set in- 
volving only A nni , C nrn , and P nm - Combining equations (34) and (36) with the defini- 
tions for & nm , 3> nm , and ®' nm reveals that B nm and D nm satisfy a set of 

homogeneous equations. Thus, if the solution for B nm and D nm is to be nontrivial, 
the determinant of the coefficients of these terms must vanish giving 
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Similarly, using equations (33), (35), and (37) along with the definitions for 
etc. , reveals that A nm , C nm , and P nrn also satisfy a set of homogeneous equations. 
For a nontrivial solution, the determinant of the coefficients of these terms must also 
vanish giving 

i;6Sv) + ;j„(v^v) n<„ + 1)^ kn(v) 

KjnfVi k Q a) cu 2 e p a 2 (l - i ^ k n (ik p r) 

Equations (38) and (39) may be appropriately called resonance equations since they 
uniquely determine the resonant frequencies or natural frequencies of the sphere after 
the plasma and sphere parameters (co p , u, v, K, and a) have been specified. Before 
presenting the solutions of these equations, it is worthwhile to consider the possible 
spatial configurations of the electromagnetic field. 

Solutions of equation (38) correspond to a field wherein B + 0, D #0, and 
A nrn = C nm = = 0. A study of equations (24) and (25) reveals that this field will 

have, in general, all three components of magnetic field (i. e. , h r , hg, and h^), but it 
can be shown by use of equations (8) and (11) that the electric field will have only 9 
and <p components. The pressure is zero since P nm = 0. Thus, solutions of equa- 
tion (38) generate a field that is transverse electric (TE) with respect to the the radial 
coordinate. Furthermore, it can be shown that solutions with n = 0 generate a null 
field and that solutions with m = 0 generate a field with only e^, h r , and hg compo- 
nents. Finally, it should be noted that equation (38) does not contain a term that is de- 
pendent on the acoustic velocity, which in turn depends on the electron temperature. 
Thus, the resonant frequencies of the TE modes are independent of the electron temper- 
ature. 

Solutions of equation (39) correspond to a field wherein A + 0, C * 0, P 0, 
and B nm = D nm = 0. It can be shown by use of equations (8), (11), (24), and (25) that 
these solutions correspond to a field that has, in general, all three components of the 
electric field (i. e. , e r , eg, and e^) but only the 9 and cp components of the magnetic 
field. These solutions are transverse magnetic (TM) with respect to the radial coordi- 
nate. The solutions for n = 0 generate a null field, as in the TE case; and solutions 
with m = 0 generate a field with only e r , eg, and h^ components. The resonant fre- 
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quencies of the TM modes are dependent on the electron temperature through the term 
kp, which appears only in the last term in equation (39). As the electron temperature 
approaches zero, k p approaches infinity, which makes the last term in equation (39) 
approach zero. Thus, the resonance equation for the TM modes in a cold plasma con- 
sists of the first two terms in equation (39). 


NUMERICAL RESULTS 

The resonant modes that correspond to the solutions of equations (38) and (39) will 
be denoted as TE^ nm and TM^, respectively. The subscript l denotes the Z** 1 
solution of the resonance equation, with l = 1 being the solution with the lowest resonant 
frequency. The subscript n denotes the order of the spherical Bessel functions, and m 
corresponds to the number of cyclic variations in the field as cp varies from 0 to 2tt. 
Since the resonance equations are independent of m, and m can take values from -n 
to n, the resonant modes will be (2n + l)-fold degenerate. For example, the TE.^ j _j, 
TEj j q, and TEj j j modes will have the same resonant frequency. 

Numerical solutions of the resonance equations were obtained by using an IBM 7090 
computer. The first case considered was the resonant frequencies of a dielectric sphere 
in a cold, collisionless plasma; that is, both the electron collision frequency v and the 
acoustic velocity u were set equal to zero. The results for this case are shown in fig- 
ure 2 in the form of curves of dimensionless resonant frequency coa/c against the di- 
mensionless plasma frequency o> p a/c. The calculations were carried out with the rela- 
tive dielectric constant of the sphere equal to 1, 9, and 100. 

The results show that, as u> a/c approaches infinity, which corresponds to the 

ir 

electron density going to infinity, the resonant frequency of each mode approaches a 
limiting value. This result is expected since a cold, collisionless plasma of infinite 
electron density acts as a conductor with infinite conductivity. Thus, the resonant fre- 
quencies of a dielectric sphere immersed in a plasma of infinite electron density are 
equal to those of a dielectric sphere coated with a perfectly conducting surface. These 
limiting resonant frequencies can be shown to be only functions of K and a. Some nu- 
merical values for these frequencies are given in table I. 

At moderate electron densities, the electromagnetic field penetrates a slight dis- 
tance into the plasma. The field, however, is still totally reflected by the plasma. 

Since the effective point of reflection is now located outside the surface r = a, the effect 
on the field is equivalent to an increase in the radius of the cavity; thus, the resonant 
frequencies are less than their values for infinite electron density. 

There are two cases to consider as the electron density becomes small. In the first 
case, the resonant frequency always remains less than the plasma frequency; thus, 
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Dimensionless resonant frequency, cja/c 



(b) Dielectric constant, K = 9. 



0 123456789 10 

Dimensionless plasma frequency, a)^tc 


(c) Dielectric constant, K = 100. 

Figure 2. - Dimensionless resonant frequency plotted against dimensionless plasma 
frequency. 
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TABLE I. - RESONANT FREQUENCIES FOR 


TABLE H. - RESONANT FREQUENCIES WHERE 


INFINITE ELECTRON DENSITY 


Mode 

Dimensionless resonant frequency, 
Vi (wa/c) 


2.74 

™12m 

3.87 

TE llm 

4.49 

™13m 

4.97 

TE 12m 

5.76 


PLASMA FREQUENCY EQUALS 
RESONANT FREQUENCY 


Mode 

Dimensionless resonant frequency, 
Vk (toa/c) 

™lnm 

0 

TE llm 

3.14 

TE 12m and ™21m 

4.49 

TE 13m and ™22m 

5.76 

TE 21m 

6.28 


solutions exist down to zero electron density. This condition only occurs for the TMj nm 
modes. Only the ™i2m’ an< * ™13m moc * es this c l ass are shown in fig- 

ure 2. It can be shown by using the asymptotic forms of the spherical Bessel functions 
for small argument along with equation (39) that the resonance curves for the TM 1 nrn 
modes are given by 


"5 = V r_JL±i_l 7 (40) 

c c [n(K + 1) + 1_ 

when coa/c « 1. This result, with K = 1, is in agreement with the quasi-static analy- 
sis of Budden (ref. 9). 

In the second case, the resonant frequency is equal to the plasma frequency for some 

electron density and exceeds the plasma frequency as the electron density is decreased 

from that value. This condition occurs for all TE modes and for the TM^ nm modes 

where l is greater than 1. It can be shown that the resonant frequencies where the 

resonant frequency is equal to the plasma frequency are only functions of K and a. 

Some numerical values for these frequencies are given in table n. 

When the resonant frequency is above the plasma frequency, energy can be radiated 

from the dielectric sphere so that a shock-excited field in the sphere will decay with 

time. The calculation of the resonant frequencies and damping factors for this case will 

be deferred until the effect of electron collisions is introduced. 

Tables I and n show that the total percentage change in the resonant frequency as 

u> varies from w„ = w to w =® is independent of K for each mode except for the 
P P P 

TM^m modes. The TM 1 r|rr| modes offer the possibility of measuring very low electron 
densities (see fig. 2, p. 15) since the resonance curves extend to a>p = 0. However, 
since the resonance curves of all TM lnm modes pass through u> p = 0, it may be diffi- 
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cult, in practice, to distinguish between these various modes for small values of Wp 
unless a technique can be devised to excite only one of these modes. It should also be 
noted that the model used for the plasma is not valid when w becomes too small. If any 
mode other than a TMj nm mode is used, the lower limit of the range of usable electron 
densities is nonzero. For example, if the TEjj m mode is used, the range of electron 

densities is restricted to where V K w a/c > 3. 14 (see table n). The lower limit of this 

ir 

range of densities can be decreased by increasing either K or a. The maximum value 
for K that can presently be attained, consistent with low electrical losses, is approxi- 
* mately 100. Increasing a increases the size of the sphere and thereby reduces the spa- 

tial resolution of the measurement. 

. A study of figure 2 reveals that the resonant frequency changes most rapidly with 

plasma frequency when the resonant and plasma frequencies are equal. Since wa/c is 
proportional to w^a/c for the TM 1 r|rn modes when u>a/c « 1 (eq. (40)), it can easily 
be shown that a 1-percent change in electron density corresponds to a 0. 5-percent 
change in resonant frequency. This result is independent of K. For the case where 
K = 1, the initial slopes of the resonance curves for modes other than the TMj nm modes 
correspond approximately to a 0. 25-percent change in resonant frequency for a 1-percent 
change in electron density. These sensitivities are reduced for higher values of relative 
dielectric constant. The high sensitivity of the sphere with K = 1 is offset by the fact 
that the resonance curves for low K approach their limiting values for infinite electron 
density more rapidly than the curves for higher K and thereby reduce the range of elec- 
tron densities that can be measured. 

It was pointed out in the section MODEL AND BOUNDARY CONDITIONS that the 
model under consideration does not take into account the sheath that forms over the sur- 
face of the sphere. This sheath can be approximated by assuming that there is a thin 
free-space layer between the dielectric and the plasma. The effect of this layer can 
easily be taken into account for the case where K = 1 since the sphere and the free- 
space layer then have the same relative dielectric constants. The net effect in this case 

is to increase the sphere radius by the thickness of the free-space layer. Since typical 

-5 

sheath thicknesses are of the order 10 meter and typical sphere radii can be of the 
-2 

order 10 meter, the sheath effectively increases the sphere radius by approximately 
0. 1 percent. This increase in sphere radius produces changes in the resonant frequen- 
cies of less than 0. 1 percent. The effect of the sheath when the relative dielectric con- 
stant of the sphere is other than 1 cannot be treated as readily. However, it is antici- 
pated that the effect will be of similar magnitude. 

It should be noted that the solutions for the case of zero electron temperature and 
collision frequency give rise to real values for k Q a and imaginary values for k g a and 
kpa. Thus, the field has an oscillatory behavior in the radial direction within the sphere 
and decays exponentially with radius outside the sphere. 
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Collisions between the electrons and other particles in the plasma are a source of 
energy loss to the field. A shock-excited field in the dielectric sphere varies with time 
as where the damping factor a accounts for the decay in the amplitude of the 

field with time. The resonant frequencies and damping factors when collisions are 
present were determined by replacing a> in the resonance equations by o> + i a and then 
finding simultaneous values for w and a that satisfied the equations. 

Numerical results for the case of a dielectric sphere immersed in a cold collisional 
plasma are shown for several modes in figures 3 and 4 in the form of curves of coa/c 



Dimensionless plasma frequency, w p a/c 

(c) Mode, TM Um . <d> Mode, TM 12nr 

Figure 3. - Dimensionless resonant frequency plotted against dimensionless plasma frequency. Dielectric constant, K*9. 
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Dimensionless damping factor, a/a) 





Dimensionless plasma frequency, 6j p a/c 


(c) Mode, TM Um . (d) Mode, TM 12m . 

Figure 4. - Dimensionless damping factor plotted against dimensionless plasma frequency. Dielectric 
constant, K = 9. 



and the dimensionless damping factor oc/oj against w p a/c for various dimensionless 
collision frequencies v / The relative dielectric constant of the sphere was set equal 
to 9. The results show that the resonant frequency wa/c is not greatly perturbed by 
the effect of collision frequencies v/ in the range 0 to 0. 1. The damping factor a/u> 
is, however, very sensitive to the value of v/ co^. Thus, for low collision rates, the 
electron density can be determined by measuring o>a/c as in the collisionless case, and 
then the collision frequency can be determined by measuring a/w. In practice, the Q 
of the cavity would normally be measured and then the damping factor would be computed 
from the equation 

Q = — 

2a 

When the collision frequency y/co p is in the range 0. 1 to 1, the resonant frequency 
u>a/c is reduced considerably from its value when v/co p - 0- Thus, for moderate 
collision frequencies, an iterative procedure must be used to find values for co p a/c 
and v/u> p that are in agreement with measured values of coa/c and a/w. If i//w p 
exceeds 1, the Q of the cavity becomes too small for this diagnostic technique to be 
useful. 

For the case of the TE llrn and TE^m modes (see figs. 2(a) and (b)), the damping 
factor, which is a strong function of v /^p at resonant frequencies much below the 
plasma frequency, becomes very large at frequencies above the plasma frequency and 
approaches a value which is independent of v /^p- This result is understandable since 
collisions represent the only energy loss mechanism at frequencies below the plasma 
frequency. At frequencies above the plasma frequency, however, energy can be lost 
from the sphere through both collisions and radiation. It is evident from figures 4(a) 
and (b) that, at frequencies above the plasma frequency, radiation is the dominant loss 
mechanism since the value of the damping factor is almost independent of collision fre- 
quency. Radiation losses are not important in the case of the TM^ nm modes since 
their resonant frequencies are always below the plasma frequency. 

Solutions for the case of a nonzero electron collision frequency give rise to complex 
values for k a = coa/c + i(aa/e), k a, and k a. Thus, the field has both oscillatory and 

O c p 

exponential behavior with radius within the sphere and, in most cases, is a traveling 
wave that decays extremely rapidly with radius external to the sphere. 

The final case to be considered is the effect of the electron temperature on the reso- 
nant frequencies and damping factors of the various modes. As discussed previously 
the resonant frequencies and damping factors of all TE modes are independent of the 
electron temperature, within the limits imposed by the model; thus, no further consider- 
ation of the TE modes is necessary. Figure 5(a) shows a plot of wa/c against co p a/c 
for the TMjj m mode for values of u/c of 0, 0.01, and 0. 1. The collision frequency 
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Dimensionless plasma frequency, co^alc 

(a) Mode, TM^ m ; dielectric constant, K = 9. (b) Mode, TM^; dielectric constant, K = 1. 

Figure 5. - Dimensionless resonant frequency against dimensionless plasma frequency. 

y/Wp was set equal to zero. As shown, the change in cua/c due to varying u/c is not 
great. The value of u>a/c only changes approximately 0. 25 percent in going from 

u/c = 0 to u/c = 0. 01, which is equivalent to changing the electron temperature from 0° 
5° 

to 3.7X10° K. 

A study of the resonance equation for the TM modes (eq. (39)) reveals that the last 
term, which is the only temperature -dependent term, contains the factor n(n + 1). 

Thus, it may be anticipated that the resonant frequencies of the modes with large n will 
have a more pronounced temperature dependence. This effect was verified from exten- 
sive numerical solutions of the resonance equation. It was also found that the tempera- 
ture effects could be increased by decreasing K. A typical result is shown in figure 5(b) 
for the mode with K = 1. The shift in resonant frequency is now approximately 

2. 7 percent in going from u/c = 0 to u/c = 0. 01. 

As discussed in the section MODEL AND BOUNDARY CONDITIONS, the validity of 
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the computed temperature effects is dependent on the validity of the rigidity boundary 
condition given by equation (14). An experimental study of the resonant frequencies of a 
dielectric sphere in a warm plasma should give a good understanding of the accuracy of 
this boundary condition. 


CONCLUDING REMARKS 

The resonant frequencies and damping factors of a dielectric sphere immersed in a 
warm collisional plasma were computed for various electron densities, electron colli- 
sion frequencies, and electron temperatures. The resonant frequencies and damping 
factors are sufficiently strong functions of the electron density and collision frequency 
so that a dielectric sphere can be useful as a plasma diagnostic tool. The variations of 
the resonant frequencies with electron temperature are also sufficiently large to be of 
diagnostic importance. However, the validity of these computed temperature effects is 
dependent on the accuracy of the rigidity boundary condition. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 14, 1967, 

129-01-07-05-22. 
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APPENDIX - SYMBOLS 


radius of sphere 

Re 

real part 

unit vectors 

r 

radial coordinate (see fig. 1) 

velocity of light 

T e 

electron temperature 

electric field intensity 

t 

time 

r, 8, and cp components of 

u 

acoustic velocity 

electric field intensity 

V 

electron fluid velocity 

magnetic field intensity 

a 

damping factor 

r, 8, and cp components of 

r 

ratio of specific heats with 

magnetic field intensity 


constant pressure to con- 
stant volume for electron 

imaginary part 


gas 

yTi 

e o 

electric permittivity of free 

relative dielectric constant 


space 

of sphere 

e p 

relative dielectric constant 

Boltzmann constant 


of plasma 

wave number for magnetic 

8 

azimuthal coordinate (see 

field in plasma 


fig. 1) 

free-space wave number, 

^0 

magnetic permeability of 

U)/c 


free space 

wave number for pressure 

V 

effective electron collision 

field in plasma 
mass of electron 


frequency for momentum 
transfer 

equilibrium electron number 

<p 

polar coordinate (see fig. 1) 

density 

perturbation electron 


magnetic field generating 
function 

pressure 

O) 

angular frequency 

electronic charge 

w p 

plasma frequency 
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